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We  rigorously  establish  the  rate  equations  for  vertical-cavity  semiconductor  optical  amplifiers, 
starting  from  a  general  energy  rate  equation.  Our  results  show  that  the  conventional  rate  equation 
used  so  far  in  the  literature  is  incorrect  because  of  an  inappropriate  calculation  of  the  mirror  losses. 
Our  calculations  include  the  effect  of  amplified  spontaneous  emission  and  can  be  used  to  describe 
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Vertical-cavity  semiconductor  optical  amplifiers 
(VCSOAs)  have  emerged  recently  as  one  important  family  of 
optoelectronic  devices.  Indeed,  conventional  in-plane  semi¬ 
conductor  optical  amplifiers  (SOAs),  unless  specially  de¬ 
signed,  show  polarization  sensitivity  and  high  coupling 
losses  to  optical  fibers.  Because  of  their  vertical  and  circular 
cavity  geometry,  VCSOAs  can  inherently  overcome  these 
problems  while  offering  the  possibility  of  parallel  processing 
and  testing.  They,  moreover,  have  the  advantage  of  a  low- 
noise  figure  compared  to  SOAs.1 

The  models  proposed  for  the  design  of  in-plane  SOAs 
(Refs.  2-4)  cannot  be  readily  extended  to  VCSOAs.  Indeed, 
because  of  their  small  gain  per  pass,  VCSOAs  require  mir¬ 
rors  with  high  reflectivity,  which  is  usually  obtained  by  tak¬ 
ing  advantage  of  additive  interferences  in  periodic  dielectric 
structures,  known  as  dielectric  Bragg  reflectors  (DBRs).  The 
amplitude  and  phase  of  these  mirrors  vary  strongly  with  the 
angle  and  wavelength.  In  addition,  because  of  the  strong 
feedback  provided  by  the  mirrors,  a  standing  optical  wave  is 
created  in  the  cavity.  This  requires  careful  positioning  of  the 
active  layers  inside  the  cavity  so  as  to  benefit  from  the  gain- 
enhancement  mechanism.5  Finally,  because  the  device  has  to 
be  biased  above  material  transparency,  amplified  spontane¬ 
ous  emission  (ASE)  is  emitted  and  cannot  be  neglected  close 
to  the  laser  threshold.  The  standard  way  of  calculating  the 
ASE  in  in-plane  SOAs  is  to  use  a  traveling  rate  equation  for 
the  field  intensity.4'6  This  is  justified  as  long  as  the  active 
region  is  longer  than  the  wavelength  of  the  cavity  modes.  For 
VCSOAs,  since  the  active  region  is  very  thin,  it  is  important 
to  consider  the  field  rather  than  the  powers  in  order  to  cor¬ 
rectly  take  into  account  interference  effects. 

Performance  predictions  for  VCSOAs  were  first  pro¬ 
posed  by  Tombling,1  based  on  the  work  of  Mukai  and 
Yamamoto3  using  semiconductor  rate  equations  for  calculat¬ 
ing  the  gain  of  the  structure.  Karlsson  and  Hoijer7  investi¬ 
gated  the  detection  and  amplification  characteristics  of 
VCSOAs  using  the  Fabry-Perot  approach  initially  derived 
by  Adams6  for  Fabry-Perot  laser  amplifiers  (FPLAs).  Re¬ 
cently,  Piprek,  Bjorlin,  and  Bowers  proposed  a  detailed 
analysis  of  VCSOAs,  based  on  the  rate  equation  approach.3 
The  rate  equation  and  the  Fabry-Perot  methods  are  known  to 
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give  significantly  different  predictions,  especially  in  the  gain 
saturation  regime.2  8  This  mechanism  occurs  as  high-signal 
power  is  injected  in  the  amplifier  or  close  to  the  laser  thresh¬ 
old:  because  of  the  corresponding  increase  of  stimulated 
emission,  the  carrier  density  in  the  active  region  is  reduced. 
This  decreases  the  material  gain  and  then  the  gain  of  the 
amplifier. 

It  is  then  of  great  interest  to  understand  the  origin  of  this 
discrepancy  and  to  rigorously  establish  the  rate  equations  for 
VCSOAs,  which  is  the  main  objective  of  this  letter.  Our 
results  show  that  the  conventional  rate  equation  used  so  far 
in  the  literature  is  incorrect  for  below  threshold  analysis  of 
the  mirror  losses.  Our  calculations  include  the  effect  of  am¬ 
plified  spontaneous  emission.  We  conclude  this  letter  demon¬ 
strating  that  our  results  can  be  readily  applied  to  resonant- 
cavity-enhanced  photodetectors  (RCEPDs),  which  can  be 
considered  as  VCSOAs  with  negative  material  gain. 

For  VCSOAs,  the  active  region  is  usually  made  of  thin 
quantum  wells  placed  at  the  antinode  of  the  electric  field  in 
the  cavity.  Provided  the  injection  is  uniform  in  the  cavity  and 
lateral  effects  are  neglected  (such  as  lateral  carrier  diffusion), 
the  material  gain  can  then  be  assumed  to  be  independent  of 
the  position.  This  simplifies  considerably  the  problem  of  cal¬ 
culating  both  the  amplified  spontaneous  emission  and  the 
gain  of  the  amplifier.  This  approximation  allows  us  to  use  a 
single  energy  rate  equation  for  the  photons,5'9  which  is  ex¬ 
pressed  as 

<?Wtot 

=  Gtot~  P tot .  (1) 

where  Wtot=hvVpNp  is  the  total  energy  stored  in  the  struc¬ 
ture,  h  v  the  energy  of  the  mode,  V p  the  volume  occupied  by 
the  photons,  and  Np  the  photon  density  in  the  active  region. 
The  total  generation  rate  can  be  expressed  as  the  sum  of 
three  contributions:  Gtot=  Gtot+Gtot+GsP-  where  Gttt  and 
Gtot  are  the  generation  and  absorption  rates  related  to  the 
gain  and  internal  losses  in  the  structure,  and  Gsp  the  power 
generated  inside  the  structure  by  the  spontaneous  emission 
and  directed  into  the  optical  mode.  With  these  definitions,  the 
modal  gain  can  be  expressed5,9  as  (g)  =  Gtp/(t;gWtot)  and  the 
internal  losses  as  «,  =  —  G U]l/(u The  confinement  factor 
T  is  by  definition  the  ratio  between  the  volume  V  occupied 
by  the  carriers  in  the  active  region,  and  the  volume  V p  oc- 
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cupied  by  the  photons.  It  is  related  to  the  modal  gain  accord¬ 
ing  to  T  =  {g)lg  =  VIVp  .  The  total  power  emitted  by  the  de¬ 
vice  is  given  by  f>tot=i3s+Pr—  Pj+Psp,  where  PR,  PT,  and 
P,  are  the  powers  reflected,  transmitted,  and  injected  in  the 
device,  respectively,  and  Psp  the  spontaneous  emission 
power  radiated  outside  the  cavity  and  in  the  optical  mode. 

Using  the  above  definitions  and  energy  rate  Eq.  (1),  one 
finds  the  general  rate  equation  for  the  photons5  dNp/dt 
=  (T g—  a  —  am)  v gNp+r /3SpRsp  where  am  is  the  mirror 
losses  defined  by  am  =  (PR  +  PT— Pj+  Psp)/(v  gh  vVpNp). 
The  spontaneous  emission  factor  /3sp  describes  the  fraction  of 
total  spontaneous  emission  that  radiates  into  the  considered 
optical  mode.  The  total  spontaneous  emission  rate  Rsp  is  re¬ 
lated  to  the  spontaneous  generation  rate  according  to  Gsp 
=  hvV /3spRsp.  In  order  to  take  into  account  the  power  emit¬ 
ted  by  the  device  only,  we  define  the  apparent  mirror  loss 
am  =  (PR  +  PT+  Psp)/(v ghvVpNp).  This  leads  to  the  rate 
equation  for  optical  amplifiers,  involving  the  amount  of 
power  injected  in  amplifier  P 7: 

dN„  P7 

—  =  (rg-a-am)vgNp+j^  +  r/3spRsp.  (2) 

This  expression  is  formally  the  same  as  the  one  generally 
used  in  the  literature.  However,  using  the  definition  of  am 
and  Eq.  (2)  at  steady  state  (dNp/dt  =  0),  we  find  that  the 
apparent  mirror  losses  must  be  expressed  by 


gain  g  can  be  considered  as  constant  in  the  active  region, 
which  is  a  reasonable  approximation  in  practice. 

For  the  sake  of  completeness,  let  us  give  the  analytical 
expressions  of  Gsignai  and  GASE  for  a  general  VCSOA  struc¬ 
ture  consisting  of  two  DBRs  surrounding  a  low-order  cavity 
of  total  thickness  Lc ,  in  which  is  placed  a  thin  active  region. 
The  DBRs  are  made  of  alternating  high  (n2)  and  low  ( n , ) 
index  layers  of  thickness  L2  and  L , ,  respectively,  with 
A.Bragg=4«i  Li  =  4n2L2.  The  complex  reflection  coefficients 
of  the  front  and  back  mirrors  are  defined  as  rl2 
=  P\  2exp(/i/'i  2),  where  f?i2=lrt2|2  =  Pi2  and  *Ai,2  is  the 
phase  of  the  mirror,  which  can  be  related  to  an  equivalent 
penetration  depth  Lpe2  of  the  electric  field  in  the  mirrors.5  An 
equivalent  cavity  length  can  then  be  defined  as  Leff=Lpen 
+  Lpen  +  Lc.  The  gain  per  pass  is  expressed  as  Gi  =  exp[(T^ 

With  these  definitions,  the  total  signal  gain  given  by 
G signai =  Gr  +  Gt  can  be  easily  calculated  according  to 

(l-R^l+R.GsKG-l) 

where  d>0  =  27rneffLeff(\0  is  the  single -pass  phase 

detuning,  neff  the  effective  refractive  index  of  the  cavity,  and 
\0  the  wavelength  of  the  signal  to  be  amplified. 

Calculation  of  gain  GASE  is  a  bit  more  complicated  and 
will  be  presented  elsewhere.  We  find: 


2m=(T  g-a,-) 


G  signal^*/ T  G  AgEG  sp 

(Gsignal-l)P7+(GASE-l)Gsp’ 


(3) 


and  depend  on  the  gain/absorption  of  the  active  region  and 
on  the  power  injected  and  generated  in  the  device.  This  re¬ 
lation  involves  two  gain  factors:  Gsignal=  GR  +  GT  is  the  total 
gain  of  the  amplifier,  with  the  gain  in  reflection  (transmis¬ 
sion)  mode  defined  by  GR  T=  P R  Tl P 7 .  The  ratio  between 
Psp  and  Gsp  is  defined  as  the  gain  of  the  amplified  spontane¬ 
ous  emission  G  ASK=  Psp/Gsp=  Psp/(h  vV(3spRsp).  This  fac¬ 
tor  indicates  how  much  ASE  power  is  extracted  outside  the 
cavity  for  the  fraction  of  spontaneous  emission  that  is  radi¬ 
ated  inside  the  cavity  and  in  the  optical  mode. 

The  photon  density  in  the  active  region  can  be  calculated 
using  Eq.  (2),  which  allows  us  to  calculate  the  stimulated 
recombination  rate  Rst=v ggNp  .  We  find 


1 


hvV  Tg  —  a 


[(Gsignal— l)R/+(GASE-l)Gsp], 


(4) 


Note  that  Rst  is  the  sum  of  two  factors  that  are  proportional 
to  the  product  between  the  inverse  of  the  amplifier  quantum 
efficiency1"  Tg  l{Tg  —  a,),  net  power  gain,  and  source  term. 
This  result  makes  obvious  sense:  for  the  input  signal,  the 
number  of  photons  generated  in  the  amplifier  is  proportional 
to  the  total  emitted  power  GsignalP7  minus  the  number  of 
injected  photons  proportional  to  P7 . 

To  completely  describe  the  properties  of  vertical-cavity 
devices,  a  rate  equation  for  the  carriers  is  needed  and  is 
generally  given  by  dN/dt  =  Ggen—Rrec,  where  Ggen  is  the 
rate  of  injected  electrons  and  R rec  is  the  rate  of  recombining 
electrons  per  unit  volume  in  the  active  region.5  The  rate 
equations  for  the  carriers  and  the  photons  can  be  used  to 

describe  any  vertical-cavity  device  provided  the  material 
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VG,2(,J>(U)(1 -/?,■)(! +  V^G,)2 

C/ACp -  , -  / -  / -  •  \0) 

2(1  +  V*  iG,)(l  +  V*2G,)(1  -  ^RiR2Gs) 

At  the  lasing  wavelength  (determined  by  the  condition  <l»(l 
=  27 r),  gains  Gsignal  and  GASE  become  infinite  when  G* 
=  1  / yjR  \  R2-  The  apparent  mirror  loss  given  by  relation  (3) 
converges  then  to  the  usual  factor 
=  ( 1/Leff)ln[l/V'R iR2]>  which  was  generally  used  for  the 
modeling  of  resonant-cavity  semiconductor  amplifiers.1'3  811 
Figure  1  displays  the  parameter  a*/2m  (GJ=l,Gsp=0)  and 
a*/ am  ( G t  =  I J1  /=())  versus  the  front  and  back  mirror  re¬ 
flectivities  as  solid  and  dotted  contour  lines,  respectively. 
These  factors  represent  the  error  that  is  made  using  a*  in¬ 
stead  of  am  in  rate  Eq.  (2)  at  transparency  (G s  =  1).  For  no 
input  signal  (Pj=  0)  and  for  highly  reflective  mirrors,  the 
apparent  mirror  losses  am  are  very  close  to  a*,  which  can 
then  be  safely  used  in  the  rate  equations  for  calculation  of  the 
amplified  spontaneous  emission.  However,  a  study  of  expres¬ 
sion  (3)  above  transparency  and  below  threshold  in  the  case 
where  the  input  signal  is  higher  than  the  spontaneous  emis¬ 
sion  ( Pi>Gsp )  shows  that  using  a*  instead  of  am  in  Eq.  (2) 
leads  to  an  overestimate  of  the  mirror  losses  by  a  factor  that 
can  be  as  high  as  4.  This  result  explains  why  the  photon 
density  calculated  with  the  conventional  rate  equation1'3'811 
leads  to  a  systematic  underestimate  of  the  photon  density 
compared  to  other  methods  based  on  Fabry-Perot  or  travel¬ 
ing  wave  approaches,  for  example.4'6'712  In  practice,  this  er¬ 
ror  leads  to  overestimating  the  saturation  power  of  VCSOAs 
by,  typically,  10  dB.  The  results  presented  in  this  letter  are 
very  general  and  can  be  applied  to  any  Fabry-Perot  laser 
amplifier  for  which  the  material  gain/absorption  can  be  con¬ 
sidered  as  constant  in  the  structure.  This  key  condition  is 

usually  satisfied  for  vertical-cavity  devices.  Note  that  the  ma- 
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Front  mirror  reflectivity  R] 

FIG.  1.  Parameters  ajJ/S,,,  (GJ=l,Gsp=0)  and  (GS=1,P;=0) 

plotted  as  solid  and  dotted  contour  lines,  respectively,  vs  the  front  and  back 
mirror  reflectivities. 

terial  gain  of  in-plane  Fabry-Perot  laser  amplifiers  cannot  be 
considered  as  constant  over  the  active  region.  Local  expres¬ 
sions  of  the  rate  equations  must  then  be  used,  as  was  pro¬ 
posed  by  Marcuse4  several  years  ago. 

It  is  interesting  to  observe  that  the  expression  of  the 
stimulated  recombination  rate  due  to  the  signal  in  Eq.  (4) 
corresponds  exactly  to  the  result  of  Adams6  who  calculated 
the  averaged  photon  density  in  in-plane  Fabry-Perot  laser 
amplifiers  assuming  a  constant  material  gain  throughout  the 
cavity  and  using  another  approach  based  on  traveling  wave 
equations  for  the  fields.  Our  expression  for  the  stimulated 
recombination  rate  due  to  the  ASE  is,  however,  different: 
because  of  the  thin  active  regions,  the  VCSOAs  must  be 
treated  considering  the  fields  rather  than  the  intensities. 

Note,  finally,  that  using  the  rate  equations  given  above,  it 
is  possible  to  simply  calculate  the  quantum  efficiency  rj  of  a 
RCEPD,  which  can  be  considered  as  a  VCSOA  with  negative 
gain.  The  quantum  efficiency  of  these  devices  is  defined  as 
the  ratio  between  the  power  absorbed  by  the  active  region, 
and  the  power  injected  in  the  device,  thus  rj=  —  Gtot/P; . 


Using  the  definitions  presented  above,  and  assuming  an  ab¬ 
sorption  per  pass  G,  =  exp[— rgLeff],  we  find  that  r;=  1 
—  G.  Using  the  relation  G  =  GR  +  GT  with  expression  (5),  we 
easily  retrieve  the  result  of  Unlii  and  Strite  [see  Eq.  (8)  of 
Ref.  13],  with  the  only  difference  that  the  standing  wave 
effect  is  automatically  taken  into  account  through  our  con¬ 
finement  factor  r,  which  includes  the  gain/absorption 
enhancement  factor. 

In  conclusion,  the  rate  equation  for  photons  has  been 
revisited  for  vertical-cavity  devices  below  threshold.  We 
have  demonstrated  that  the  mirror  loss  expression,  which  is 
traditionally  used  in  rate  equations  for  photons  in  resonant- 
cavity  semiconductor  optical  amplifiers,  is  incorrect  below 
threshold  because  it  does  not  satisfy  the  energy  conservation. 
This  error  leads  then  to  an  underestimate  of  the  photon  den¬ 
sity  in  the  cavity,  hence,  to  an  underestimate  of  the  stimu¬ 
lated  recombination  rate,  and  then  to  an  overestimate  of  gain 
saturation  effects  for  these  amplifiers.  The  correct  expression 
for  the  mirror  losses  has  been  derived  and  includes  the  am¬ 
plified  spontaneous  emission.  Finally,  we  have  presented  a 
unified  way  of  describing  vertical-cavity  devices  by  photon 
rate  equations,  which  can  be  applied  to  any  device  with 
stimulated  gain  or  absorption  like  VCSOAs  or  RCEPDs,  the 
only  condition  being  that  the  material  properties  of  the  gain/ 
absorbing  medium  are  constant  throughout  the  active  region, 
which  is  generally  reasonable  for  this  kind  of  device. 

This  work  was  supported  by  DARPA  through  the  Army 
Research  Office. 
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